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ABSTRACT 

We present various infrared two-color diagrams (2CDs) for AGB stars, post-AGB stars, and Plane¬ 
tary Nebulae (PNe) and investigate possible evolutionary tracks. We use catalogs from the available 
literature for the sample of 4903 AGB stars (3373 O-rich; 1168 C-rich; 362 S-type), 660 post-AGB 
stars (326 post-AGB; 334 pre-PNe), and 1510 PNe in our Galaxy. For each object in the catalog, 
we cross-identify the IRAS , AKARI , MSX , and 2MASS counterparts. The IR 2CDs can provide 
useful information about the structure and evolution of the dust envelopes as well as the central stars. 

To find possible evolutionary tracks from AGB stars to PNe on the 2CDs, we investigate spectral 
evolution of post-AGB stars by making simple but reasonable assumptions on the evolution of the 
central star and dust shell. We perform radiative transfer model calculations for the detached dust 
shells around evolving central stars in the post-AGB phase. We find that the theoretical dust shell 
model tracks using dust opacity functions of amorphous silicate and amorphous carbon roughly coin¬ 
cide with the densely populated observed points of AGB stars, post-AGB stars, and PNe on various 
IR 2CDs. Even though some discrepancies are inevitable, the end points of the theoretical post-AGB 
model tracks are generally converged to the region of the observed points of PNe on most 2CDs. 

Subject headings: stars: AGB and post-AGB - circumstellar matter - infrared: stars - dust, extinction 
- radiative transfer. 


1. INTRODUCTION 

Asymptotic giant branch (AGB) stars are generally 
classified as O-rich (M-type) or C-rich (C-type) based on 
the chemistry of the photosphere and/or the outer en¬ 
velope. When an AGB star of intermediate mass range 
goes through C dredge-up process and thus the C/O ra¬ 
tio is larger than 1, O-rich dust grain formation ceases 
and the star may become a carbon star (e.g., Iben 1981; 
Chan & Kwok 1990). S stars are generally regarded as 
intermediate between M-type and carbon stars in their 
properties (Lloyd Evans & Little-Marenin 1999). How¬ 
ever, this M-S-C evolutionary sequence can be different 
depending on the mass and metallicity (e.g., Gronewegen 
et al. 1995; Suh 2014). Dust envelopes around AGB stars 
are believed to be a main source of interstellar dust. The 
outflowing envelopes around AGB stars are very suitable 
places for massive dust formation. Nearly all AGB stars 
can be identified as Long-Period Variables (LPVs). The 
AGB phase of the LPV is characterized by dusty stel¬ 
lar winds with high mass-loss rates (10 -8 — 10 _4 MQ/j/r; 
e.g., Loup et al. 1993). 

As the star leaves the AGB phase, its mass-loss rate de¬ 
creases significantly and the star may become hot enough 
to ionize its circumstellar material. When mass-loss re¬ 
duces the mass of the remaining H-rich envelope below ~ 
10 _2 Mq (the exact value depends on the initial mass), 
the stellar envelope begins to shrink and the effective 
temperature starts to increase (e.g., Schonberner 1983) 
until the central star is hot enough (~ 30000 K) to ionize 
the circumstellar nebula. If the temperature increases on 
a timescale shorter than the dispersion time of the mat¬ 
ter previously ejected by the star, a planetary nebula 
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(PN) will appear. Planetary nebulae (PNe) are believed 
to be a common end-point of stellar evolution for a large 
fraction of all stars between 1 and 10 Mg. 

The intermediate phase between the end of the AGB 
phase and the PN phase is called the post-AGB phase. 
During the the post-AGB phase, the dust shell formed 
in the AGB phase detaches from the central star and 
becomes optically thin after a few hundred years (e.g., 
Hrivnak et al. 1989). There could be some post-AGB 
objects that have left the AGB stage but will evolve to 
the white dwarf stage without ever becoming a PN (e.g., 
Zuckernran 1978; Szczerba et al. 2007). The evolution¬ 
ary lifetime of a PN is critically dependent on the core¬ 
mass of its central star (e.g., Renzini 1981). A low-mass 
star evolves too slowly to become a PN and a high-mass 
star is luminous too briefly to be detected. Therefore, 
only stars with core-masses in a narrow mass range (~ 
0.6 M 0 ) would be seen as PNe (e.g., Kwok 2000). 

The Infrared Astronomical Satellite (IRAS) Point 
Source Catalog (PSC) (version 2.1) provided useful pho¬ 
tometric data in four bands (12, 25, 60, and 100 /mi). 
In characterizing the circumstellar environments of AGB 
and post-AGB stars, the two-color diagram (2CD) in the 
IRAS PSC has been useful (e.g., Suh & Kwon 2011). 
The IRAS Low Resolution Spectrograph (LRS; A = 
8—22 /im) data are useful to identify important features 
of dust grains. Kwok et al. (1997) used IRAS LRS to 
identify the class E (the 10 /tin silicate feature in emis¬ 
sion), class A (the 10 /im silicate feature in absorption), 
class C (the 11 /im SiC dust emission), and class P (the 
11.3 /tin or 12.5 /mr emission features that are attributed 
to polycyclic aromatic hydrocarbon (PAH)). 

The Midcourse Space Experiment (MSX; Egan et al. 
2003) surveyed the Galactic plane in four bands (8.28, 
12.13, 14.65, and 21.34 /im). AKARI (Murakami et al. 
2007) provided PSC data in two bands (9 and 18 /an) 
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from an all-sky survey. The PSC from the two micron all 
sky survey (2MASS; Cutri et al. 2003) contains accurate 
positions and fluxes in J (1.25 /xm), H (1.65 /xm), and 
I\ s (2.17 /xm) bands. 

In this paper, we present various IR 2CDs for AGB 
stars, post-AGB stars, and PNe and investigate possi¬ 
ble evolutionary tracks. We use catalogs from the avail¬ 
able literature for the sample stars in our Galaxy. For 
each object in the catalog, we cross-identify the IRAS, 
AKARI , MSX, and 2MASS counterparts. To investi¬ 
gate possible evolutionary tracks from AGB stars to PNe 
on the 2CDs, we make simple but reasonable theoretical 
models on the evolution of the central star and dust shell 
in the post-AGB phase. We perform radiative transfer 
model calculations for detached dust shells around evolv¬ 
ing central stars in the post-AGB phase. We compare the 
theoretical model tracks with observations of AGB stars, 
post-AGB stars, and PNe and discuss possible evolution¬ 
ary tracks on the IR 2CDs. 

2. SAMPLE STARS AND IR TWO-COLOR DIAGRAMS 

We use catalogs from the available literature for the 
sample of 4903 AGB stars (3373 O-AGB; 1168 C-AGB; 
362 S-type), 660 post-AGB stars (326 post-AGB; 334 
pre-PN), and 1510 PNe. For each object in the catalog, 
we cross-identify the IRAS PSC, AKARI PSC, MSX 
PSC, and 2MASS counterparts as we will explain in sub¬ 
sections. 

In Table 1, we list the reference, total number of ob¬ 
jects, and numbers of the cross-identified IRAS PSC, 
AKARI PSC, MSX PSC, and 2MASS counterparts. 
We also list the number of sources according to the IRAS 
LRS classification. Class E and A objects show O-rich 
dust features and class C and P objects show C-rich dust 
features (see section 1). 

There are considerable differences in the angular reso¬ 
lutions of different data sets. The angular resolutions for 
the IRAS, MSX, AKARI, and 2MASS are 0.75x4.5- 
4.6', 18.3", 9.4" (in 9 /xin band), and 2", respectively. 
The differing spatial resolutions may have impacts on 
the photometry used in the 2CDs. The 2MASS flux pre¬ 
sumably includes only the central star, while the IRAS 
fluxes may include the extended outer envelope. 

Because IRAS has very low angular resolution, the 
survey regions may suffer from confusion problems. So 
the 2MASS counterpart obtained from the IRAS PSC 
position can be a different object. The 2MASS counter¬ 
part obtained from the position data from AKARI which 
has a higher resolution would be much more reliable. 

2.1. AGB stars 

Suh & Kwon (2011) presented a catalog of AGB stars 
for 3003 O-rich, 1168 C-ricli, and 362 S-type objects in 
our Galaxy. For the catalog, they compiled previous 
works with verifying processes from the sources listed in 
the IRAS PSC. Kwon & Suh (2012) presented a revised 
sample of 3373 0-ricli AGB stars. 

By using the position information in version 2.1 of the 
IRAS PSC, we cross-identify the AKARI and MSX 
counterparts by finding the nearest source within 30" for 
each object. For 2MASS, we find the closest counterpart 
in the position within 30" using the position information 
of the cross-identified AKARI PSC source. 


2.2. Post-AGB stars 

The intermediate phase between the end of the AGB 
phase and the PN phase is called the post-AGB phase 
(e.g., Szczerba et al. 2007). It was formerly also referred 
to as the proto-PN phase or the pre-PN phase (e.g., Zuck- 
errnan 1978; Hrivnak et al. 1989; Kwok 2000; Kohoutek 
2001 ). 

For this paper, we use the list of 326 ‘very likely post- 
AGB objects’ from the catalog of post-AGB stars by 
Szczerba et al. (2007) and 334 pre-PNe from the cat¬ 
alogue of Galactic Planetary Nebulae (Updated Version 
2000; Kohoutek 2001). We regard the 660 objects in both 
catalogs (post-AGB objects and pre-PNe) as post-AGB 
stars. 

Because AGB stars and post-AGB stars share some 
similarities, they could have been duplicately classified. 
Table 2 lists them. For S type AGB stars, there are no 
duplicated objects. 

For post-AGB stars, we use the IRAS PSC, MSX, 
and 2MASS counterparts listed in Szczerba et al. (2007). 
We cross-identify the AKARI counterparts by finding 
the nearest source within 30" using the position infor¬ 
mation in Szczerba et al. (2007) for each object. 

For pre-PNe, we use the IRAS PSC counterparts listed 
in Kohoutek (2001). We cross-identify the AKARI and 
MSX counterparts by finding the nearest source within 
30" using the position information in Kohoutek (2001). 
For the 2MASS counterpart, we find the closest source 
in the position within 30" using the position information 
of the cross-identified AKARI PSC source. 

2.3. Planetary Nebulae 

For this paper, we use the catalogue of Galactic Plan¬ 
etary Nebulae (Updated Version 2000; Kohoutek 2001). 
This catalogue contains 1510 objects classified as galactic 
PNe up to the end of 1999. Valuable information about 
the objects including their identification charts and the 
lists of references can be found in the Strasbourg-ESO 
Catalogue of Galactic Planetary Nebulae (Acker et al. 
1992; SECGPN) and in its First Supplement (Acker et 
al. 1996). Kerber et al. (2003) presented more accurate 
position data for 1312 PNe. 

There have been significant new discoveries of PNe af¬ 
ter Kohoutek (2001). The most important contribution 
would be the Macquarie/AAO/Strasbourg Ha (MASH) 
catalogue (Parker et al. 2006) which represents the 
largest incremental increase in Galactic PN numbers to 
date. However, MASH PNe are typically more evolved, 
obscured, and of lower surface brightness than those 
found in most previous surveys. Therefore, we do not 
use the newly discovered PNe because they are not very 
useful for the purpose of this work. 

For PNe, we use IRAS PSC counterpart listed in Ko¬ 
houtek (2001). We cross-identify the AKARI and MSX 
counterparts by finding the nearest source within 30" 
using the position information in Kohoutek (2001) and 
Kerber et al. (2003). For the 2MASS counterpart, we 
find the closest source in the position within 30" using 
the position information of the cross-identified AKARI 
PSC source. 

2.4. IR Two-Color Diagrams 
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TABLE 1 

Sample of AGB stars, post-AGB stars, and PNe 


Class 

Reference 

Number 

IRAS PSC 

AKRAI PSC 

MSX PSC 

2MASS 

A* 

E* 

C* 

p* 

O-AGB 

Kwon Sz Suh (2012) 

3373 

3373 

2708 

1753 

2708 

186 

1125 

0 

44 

C-AGB 

Suh & Kwon (2011) 

1168 

1168 

1012 

687 

1012 

0 

0 

713 

11 

S-AGB 

Suh Sz Kwon (2011) 

362 

362 

336 

128 

336 

0 

29 

1 

0 

post-AGB 

Szczerba et al. (2007) 

326 

236 

254 

117 

318 

4 

14 

1 

3 

Pre-PNe 

Kohoutek (2001) 

334 

326 

253 

177 

253 

15 

24 

6 

4 

PNe 

Kohoutek (2001) 

1510 

927 

808 

414 

808 

1 

15 

0 

13 


*: The IRAS LRS class. 


TABLE 2 

Duplicate classification of AGB and post-AGB stars 



O-AGB 

C-AGB 

post-AGB 

Pre-PNe 

O-AGB 

3373 

0 

34 

74 

C-AGB 

0 

1168 

12 

23 

post-AGB 

34 

12 

326 

100 

Pre-PNe 

74 

24 

100 

334 


Infrared fluxes from the IRAS (PSC; 12, 25, 60, and 
100 pm), MSX (PSC; 8.28, 12.13, 14.65, and 21.34 pm), 
AKARI (PSC; 9 and 18 /mi) and 2MASS (1.25, 1.65, 
and 2.17 /mi) are available for the sample of AGB stars, 
post-AGB stars, and PNe as discussed in previous sub¬ 
sections. The large number of observations can be used 
to form various 2CDs that can be compared with theo¬ 
retical model predictions. 

The color index is defined by 

71 r 7i f oc1 F\ 2 /ZMC\2 

m \i - = 2.5 log 10 Fa ! ZMCa (1) 

where ZMCxi means the zero magnitude calibration at 
given wavelength (Az) (see Suh & Kwon 2011 for details). 

Figures 1, 2, 3, and 4 show IRAS, NIR, AKARI, and 
MSX 2CDs for AGB stars, post-AGB stars, and PNe 
compared with theoretical models. The small symbols 
are observational data and the lines with symbols are 
theoretical model calculations. For each figure, the up¬ 
per panel plots all observed objects and theoretical model 
tracks for AGB stars. The lower panel plots observed ob¬ 
jects of post-AGB stars and PNe and theoretical model 
tracks for post AGB stars. We will discuss the theoret¬ 
ical models in sections 3 and 4. We will compare the 
theoretical models with the observations on the IR 2CDs 
in section 5. 

For AGB stars, we plot only those objects with good 
quality (q = 3) at any wavelength on all 2CDs using the 
IRAS, AKARI, 2MASS, and MSX data. On the other 
hand, we plot all objects for post-AGB stars and PNe. 
For the post-AGB stars and PNe with good quality (q = 
3) data, we plot them by thicker symbols. For the post- 
AGB stars and PNe which have the IRAS LRS data, we 
indicate the IRAS LRS class by corresponding symbols. 
Class E and A objects show O-rich dust features and 
class C and P objects show C-rich dust features. The 42 
post-AGB stars which are known to be binary systems 
(Szczerba et al. 2007) are marked by the green hexagonal 
symbol. 

For AGB stars, the objects in the upper-right regions 
on any 2CDs have thicker dust shells with large optical 
depths. Generally, the [12] —[25], [14.65] —[21.34], and 
[25] —[60] colors for PNe are redder than those for AGB 


stars with thick dust shells. The [8.28] —[14.65] colors for 
PNe are comparable to those for AGB stars with thick 
dust shells. However, K s — [12] and [9] —[12] colors for 
PNe are bluer than those for AGB stars with thick dust 
shells. These effects are clearly displayed in the 2CDs. 

The locations of PNe look to be converged toward upper- 
right regions on the two 2CDs using [25] —[60] versus 
[12] —[25] and [8.28]-[14.65] versus [14.65] — [21.34]. But 
for the 2CD using [12]—[25] versus K s — [ 12], the conver¬ 
gence occurs toward the upper-middle region. For the 
2CD using [9] —[12] versus [12] — [25], the convergence oc¬ 
curs toward the right region. 

3. THEORETICAL DUST SHELL MODELS 

To investigate the spectral evolution of AGB and post- 
AGB stars on the IR 2CDs, we use radiative transfer 
models for spherically symmetric dust shells. We assume 
that the central star emits blackbody radiation (see sec¬ 
tion 6.1 for a discussion on the assumptions). 

3.1. Radiative transfer model calculations 

In this work, we use the ra¬ 
diative transfer code RADMC-3D 

(http://www.ita.uni-heidelberg.de/~dullemond/software/radmc-3d 
to calculate the model SEDs for dust shells around AGB 
stars and pos-AGB stars. RADMC-3D is based on 
the Monte Carlo simulation method of Bjorkman & 

Wood (2001) for investigating dust continuum radiative 
transfer processes. For a spherically symmetric dust 
shell, we have compared the model results obtained 
by the RADMC-3D code with those obtained by the 
radiative transfer codes CSDUST3 (Egan et al. 1988) 
and DUSTY (Ivezic & Elitzur (1997). All of these three 
codes produced essentially identical results for the same 
model parameters. 

We assume a spherically symmetric dust shell around 
a single star. We use similar schemes for dust density 
distribution as those used by Suh & Kwon (2013). We 
assume that the dust density distribution is continuous 
from the inner radius (Ri n ) to the outer radius ( R ou t ) for 
a spherically symmetric dust shell. For the dust density 
distribution, we use the simple power law equation 

p(r) = Pi n (r/R in )~ 2 , (2) 

where pi n is the dust density at Ri n of the dust shell. 

The mass-loss rate is given by 

M — 47 TV pV eX p — 1 77 R~ n Piri.'Vexp , ( 3 ) 

where p m is the material density at Ri n which is given 

by 


Pmd — Pin: 


( 4 ) 
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Fig. 1.— IRAS 2CDs for AGB stars, post-AGB stars, and planetary nebulae compared with theoretical models. See sections 3-5 for 
discussion on the theoretical models. 


where S is the dust-to-gas ratio which is assumed to be 
0.01 (e.g., Suh 2014). The mass-loss rate (M) and dust 
shell expansion velocity ( v exp ) remain constant. 

The inner shell dust temperature ( T c ) is the dust tem¬ 
perature at Ri n . The model SEDs are sensitively depen¬ 
dent on Ri n (or T c ). For all AGB and post-AGB stars, we 
assume that the dust formation (or condensation) tem¬ 
perature (T c ) is 1000 K (e.g., Suh 1999, 2004). 

If dust formation ceases when the post-AGB phase 
starts, the dust shell begins to detach (i.e., Ri n increases 
for colder T c ). For all of the models, we use the outer 
radius ( R ou t ) of the dust shell so that the dust temper¬ 


ature at Rout is lower than 30 K which is comparable to 
the temperature of interstellar medium. 

We assume that all dust grains are spherical with a 
uniform radius of 0.1 /im and the scattering is assumed to 
be isotropic. We choose 10 //m as the fiducial wavelength 
that sets the scale of the dust optical depth (r 10 ). 

For the central star, we assume that it emits blackbody 
radiation for a given temperature and luminosity. 

3.2. AGB stars 

For AGB stars, we use similar assumptions for the dust 
shell models as those used for Suh & Kwon (2011). For 
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Fig. 2.— NIR -IRAS 2CDs for AGB stars, post-AGB stars, and planetary nebulae compared with theoretical models. 


O-rich stars, we use the optical constants of warm and 
cold silicate grains derived by Suli (1999). We compute 
models for eleven dust optical depths (ti 0 = 0.005, 0.01, 
0.05, 0.1, 0.5, 1, 3, 7, 15, 30, and 40). We use the warm 
silicate dust grains for low mass-loss rate O-rich AGB 
stars (7 models with t w < 3) and the cold silicate grains 
for high mass-loss rate O-rich AGB stars (4 models with 
no > 3). For the central star, we assume that the lumi¬ 
nosity is 10 4 Lq and the stellar blackbody temperature 
is 2500 K for no < 3 and 2000 K for no > 3. 

For C-rich stars, we use the optical constants of amor¬ 
phous carbon (AMC) grains derived by Suh (2000). We 
perform the model calculations for seven dust optical 
depths (no = 0.01, 0.1, 1, 2, 3, 5, and 7). For the central 
star, we assume that the luminosity is 10 4 Lq and the 
stellar blackbody temperature is 2300 K for no < 0.1 


and 2000 K for no > 0.1. 

In Figures 1-4, the theoretical model tracks for AGB 
stars are plotted on IR 2CDs. Suh & Kwon (2011) pre¬ 
sented more various model tracks for AGB stars using 
more dust species. 

4. SPECTRAL EVOLUTION OF POST-AGB STARS 

To understand possible evolutionary tracks from AGB 
stars to PNe on 2CDs, we investigate the spectral evolu¬ 
tion of post-AGB stars using the radiative transfer model 
as described in section 3.1. We make simple but reason¬ 
able assumptions on the evolution of properties of the 
dust shell and the central star to find theoretical evolu¬ 
tionary tracks of post-AGB stars on 2CDs. 

The SEDs of typical post-AGB stars show two distinct 
components: the cooler component corresponds to the 
remnant of the AGB dust shell, and the warmer compo- 
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Fig. 3.— AKARI-IRAS 2CDs for AGB stars, post-AGB stars, and planetary nebulae compared with theoretical models. 


nent corresponds to the photosphere of the central star 
(e.g., Volk & Kwok 1989). The clear separation of the 
two flux components would be due to the detached dust 
shell because the AGB mass-loss terminated some time 
ago. The dynamical time scale (after the cessation of 
AGB mass-loss) can be estimated by radiative transfer 
models. The typical time scales of the post-AGB phase 
range from several hundred to a thousand years (e.g., 
Hrivnak et al. 1989). 

4.1. The central star 

For the central star, we assume a fixed stellar luminos¬ 
ity and changing stellar blackbody temperatures during 
the post-AGB phase. We use relevant model parameters 
from the theoretical models of post-AGB stars presented 
by van Hoof et al. (1997). Their models are based on the 
works by Wood (1990), Vassiliadis & Wood (1993, 1994), 


Blocker & Schonberner (1991), and Blocker (1995a,b). 
Many assumptions have to be made about the mass- 
loss rate, the exact time of the end of the AGB phase, 
and the outflow velocity. The duration of the post-AGB 
phase may be as long as 10 4 -10 5 years for the less mas¬ 
sive objects, but more massive objects may last only a 
few decades or centuries in the post-AGB phase (Blocker 
1995b). 

In this work, we assume that the AGB-type mass-loss 
stops (i.e., the AGB phase ends) when the central star 
has reached an inferred pulsation period of P a = 100 
days (Blocker 1995b). When the star has subsequently 
reached an inferred pulsation period of Pb = 50 days 
the post-AGB phase (mass-loss) starts. The mass-loss 
rate changes smoothly between the AGB and post-AGB 
phases. 
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PNe (1510; Kohoutek 2001) 
pre-PNe (334; Kohoutek 2001) 
Post-AGB (326; Szczerba 2007) 


O-rich AGB stars (3373) 
C-rich AGB stars (1168) 
S stars (362) 



[14.65]-[21.34] 



Fig. 4.— MSX 2CDs for AGB stars, post-AGB stars, and planetary nebulae compared with theoretical models. 


For P a = 100 days and P& = 50 days, the calculated 
evolutionary time scales of the post-AGB phase show 
wide variations (about 38 - 8700 years; see van Hoof et 
al. 1997) depending on the core-mass ( M core ) and post- 
AGB mass-loss rates. If we assume earlier ends of the 
AGB phase (a larger value of P a ), the time scale can 
be significantly increased. Comparison of the prelimi¬ 
nary dust shell model (see section 4.2) results with the 
observations on the IR 2CDs indicated that the theoret¬ 
ical post-AGB phase needs to be rather short (within a 
thousand years) to match the observations. 

We restrict ourselves to two different core-mass models 
for this paper. For low-mass O-rich stars and C-rich 
stars, we use M core = 0.605 M 0 and L* = 6310 Lq 
with the post-AGB mass-loss rate at 5 times the standard 
value. For high-mass O-rich stars, we use M core = 0.696 


Mq and L* = 11610 Lq with the post-AGB mass-loss 
rate at the standard value. The corresponding model 
parameters are adopted from Table 1 in van Hoof et al. 
(1997) for the evolution of the central star in the post- 
AGB phase. The time scales of the post-AGB phase are 
832 and 174 years for the low core-mass and high core¬ 
mass models, respectively. 

For the low core-mass model, Table 3 presents five base 
models for the dust shell. For the high core-mass model, 
Table 4 presents three base models for the dust shell. 
We will discuss the dust shell models in section 4.2. For 
each dust shell model, the time evolution of the surface 
temperature (T e //) of the central star at eight points in 
the post-AGB phase as well as the fixed stellar luminosity 
(L*) are presented. 

4.2. The dust shell 



















TABLE 3 

Theoretical model parameters of the central star and dust shell for low core-mass post-AGB stars. 


M core = 0.605 Mq, L* = 

6310 Lq 

, Dust: Silicate, 

M=(EM:3.73 X 10 -9 ; ML:1.21 X 10" 

~ 7 : M:1.27 X 10“ 5 ) M 0 /yr 

Base Model 

Teff (K) 

tpa (years) Ri n (AU) 

Pin, (g/cm 3 ) 

T10 

Tc (K) 

Derived Model(s) 

EM 

6042 

0 

1.70E+01 

1.93E-21 

1.00E-03 

1000 

EM2(75): v exp =75 

Vexp = 15 

6500 

113 

3.75E+02 

3.98E-24 

4.59E-05 

205 

EM2(150): Vexp—150 

Ti=0.001 

7000 

178 

5.80E+02 

1.66E-24 

2.96E-05 

174 

Ti=0.001 


10000 

241 

7.80E+02 

9.18E-25 

2.20E-05 

168 



15000 

325 

1.05E+03 

5.10E-25 

1.64E-05 

164 



20000 

452 

1.45E+03 

2.66E-25 

1.18E-05 

152 



25000 

625 

1.99E+03 

1.40E-25 

8.53E-06 

138 



30000 

832 

2.65E+03 

7.95E-26 

6.44E-06 

126 


ML 

6042 

0 

1.90E+01 

5.00E-20 

2.90E-02 

1000 

ML2(75): v exp =75 

Vexp = 15 

6500 

113 

3.77E+02 

1.27E-22 

1.47E-03 

205 

ML2(150): v e xp=150 

Tj—0.029 

7000 

178 

5.82E+02 

5.32E-23 

9.51E-04 

174 

ML2(300): v exp =300 


10000 

241 

7.82E+02 

2.95E-23 

7.09E-04 

168 

Ti =0.029 


15000 

325 

1.05E+03 

1.65E-23 

5.32E-04 

164 



20000 

452 

1.45E+03 

8.59E-24 

3.82E-04 

152 



25000 

625 

2.00E+03 

4.53E-24 

2.78E-04 

138 



30000 

832 

2.65E+03 

2.57E-24 

2.08E-04 

126 


M 

6042 

0 

2.47E+01 

3.10E-18 

2.35E+00 

1000 

M2(30): Vexp—30 

Vexp = 15 

6500 

113 

3.82E+02 

1.29E-20 

1.51E-01 

241 

n= 2.35 

n=2.35 

7000 

178 

5.89E+02 

5.46E-21 

9.88E-02 

198 



10000 

241 

7.88E+02 

3.04E-21 

7.36E-02 

190 



15000 

325 

1.05E+03 

1.70E-21 

5.48E-02 

180 



20000 

452 

1.46E+03 

8.91E-22 

3.99E-02 

162 



25000 

625 

2.00E+03 

4.70E-22 

2.88E-02 

144 



30000 

832 

2.66E+03 

2.67E-22 

2.17E-02 

130 


Mcore = 0.605 Mq, L* = 

6310 Lq 

, Dust: Amorphous Carbon, M— 

= (CL:7.85 x 10 _s ; C:1.68 X 10“ 5 ) M Q /yr 

Base Model 

Teff (K) 

tpa (years) Ri u (AU) 

Pin (g/cm 3 ) 

T10 

Tc (K) 

Derived Model(s) 

CL 

6042 

0 

1.78E+01 

3.70E-20 

1.00E-02 

1000 

CL2(30): Vexp=30 

Vexp = 15 

6500 

113 

3.75E+02 

8.32E-23 

4.77E-04 

298 

CL2(150): Vexp—150 

T^—0.01 

7000 

178 

5.81E+02 

3.47E-23 

3.08E-04 

255 

CL2(300): v e xp=300 


10000 

241 

7.80E+02 

1.93E-23 

2.30E-04 

241 

T;=0.01 


15000 

325 

1.05E+03 

1.07E-23 

1.72E-04 

222 



20000 

452 

1.45E+03 

5.59E-24 

1.24E-04 

197 



25000 

625 

2.00E+03 

2.94E-24 

8.99E-05 

175 



30000 

832 

2.65E+03 

1.67E-24 

6.77E-05 

156 


c 

6042 

0 

1.83E+01 

7.50E-18 

2.11E+00 

1000 

C2(30): Vexp—30 

Vexp = 15 

6500 

113 

3.76E+02 

1.78E-20 

1.03E-01 

311 

T; = 2.11 

Ti = 2.11 

7000 

178 

5.82E+02 

7.43E-21 

6.67E-02 

266 



10000 

241 

7.81E+02 

4.12E-21 

4.96E-02 

251 



15000 

325 

1.05E+03 

2.29E-21 

3.71E-02 

232 



20000 

452 

1.45E+03 

1.20E-21 

2.68E-02 

207 



25000 

625 

2.00E+03 

6.31E-22 

1.94E-02 

182 



30000 

832 

2.65E+03 

3.57E-22 

1.45E-02 

162 



Dust shell parameters are presented for base models (see section 4.2). The unit for v exp is km/sec. Note that the model tracks for 
EM2(75), EM2(150), ML2(150) ,ML2(150), ML2(300), and CL2(300) are plotted only on the IRAS 2CD (Figure 1). 


For the dust shell, we use dust radiative transfer model 
described in section 3.1. We use the same dust opacity 
functions as those used for AGB stars, we use the opti¬ 
cal constants of cold silicate dust grains derived by Suh 
(1999) for O-rich post-AGB stars. For C-ricli post-AGB 
stars, we use the optical constants of amorphous carbon 
(AMC) dust grains derived by Suh (2000). 

For C-rich post-AGB stars, we use AMC rather than 
graphite which is widely used for interstellar dust or 
young stellar objects. It is believed that AMC is the 
main dust component for C-rich post-AGB stars (e.g., 
Hony et al. 2003; Cerrigone et al. 2009) as well as for 
C-rich AGB stars. Unlike AGB stars, C-rich post-AGB 
stars show dust features from hydrogenated amorphous 
carbon and PAH as well as SiC and MgS (e.g., Justtanont 
et al. 1996a). 


The typical dust shell expansion velocity ( v exp ) is 
about 15 km/sec for AGB stars (e.g., Loup et al. 1993; 
Suh 2014). Even though it is not clear whether v exp can 
be increased during the AGB phase or post-AGB phase, 
v exp could be higher in the post-AGB phase. Slijkhuis 
et al. (1991) measured the outflow velocities up to v exp 
= 400 km/sec for a post-AGB star (IRAS 08005-2356). 
We use various possible values of v exp (15, 30, 75, 150, 
and 300 km/sec) for post-AGB stars. 

Uncertainty on the mass-loss rate remains. A thermal 
pulse in the AGB phase can change the mass-loss rate 
abruptly (e.g., Renzini 1981; Wood 1990; Vassiliadis & 
Wood 1993, 1994) from the continuous analytical laws 
which were used for the central star (see section 4.1). 
The effect of the thermal pulse on T e ff would be mi¬ 
nor, but the effect on L* could be significant. However, 
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TABLE 4 

Theoretical model parameters of the central star and dust shell for high core-mass post-AGB stars. 


Mcore — 

0.696 M 0 , L* = 

11610 L 0 

, Dust: Silicate, M=(OL:8.64 x 10 5 ; 0:3.40 x 10 

' 4 ; T:8.38 X 10" 4 ) M Q /yr 

Base Model T e ff (K) 

tp a (years) R in (AU) 

Pin (g/cm 3 ) 

TIO 

T c (K) 

Derived Model 

OL 

6846 

0 

2.54E+01 

2.00E-17 

1.55E+01 

1000 

OL2(75): v exp =75 

Vexp = 15 

7000 

29 

1.17E+02 

9.40E-19 

3.38E+00 

558 

Ti= 15.5 

n =15.5 

7500 

62 

2.22E+02 

2.63E-19 

1.79E+00 

415 



10000 

84 

2.91E+02 

1.52E-19 

1.36E+00 

392 



15000 

111 

3.77E+02 

9.10E-20 

1.05E+00 

369 



20000 

129 

4.34E+02 

6.86E-20 

9.13E-01 

354 



25000 

151 

5.03E+02 

5.10E-20 

7.86E-01 

334 



30000 

174 

5.76E+02 

3.89E-20 

6.86E-01 

314 


O 

6846 

0 

3.00E+01 

5.65E-17 

5.20E+01 

1000 

02(75): v exp =75 


7000 

29 

1.22E+02 

3.43E-18 

1.21E+01 

501 

75=52 

Ti =52 

7500 

62 

2.26E+02 

9.94E-19 

6.48E+00 

373 



10000 

84 

2.96E+02 

5.81E-19 

4.96E+00 

343 



15000 

111 

3.81E+02 

3.50E-19 

3.84E+00 

321 



20000 

129 

4.38E+02 

2.65E-19 

3.34E+00 

310 



25000 

151 

5.08E+02 

1.97E-19 

2.89E+00 

296 



30000 

174 

5.81E+02 

1.51E-19 

2.53E+00 

283 


T 

6846 

0 

3.75E+01 

8.90E-17 

1.02E+02 

1000 

T2(75): v exp =75 

Vexp = 15 

7000 

29 

1.29E+02 

7.49E-18 

2.97E+01 

485 

Ti =102 

71=102 

7500 

62 

2.34E+02 

2.29E-18 

1.65E+01 

360 



10000 

84 

3.03E+02 

1.36E-18 

1.26E+01 

320 



15000 

111 

3.89E+02 

8.28E-19 

9.87E+00 

290 



20000 

129 

4.46E+02 

6.30E-19 

8.60E+00 

276 



25000 

151 

5.15E+02 

4.71E-19 

7.42E+00 

263 



30000 

174 

5.88E+02 

3.62E-19 

6.50E+00 

252 



Dust shell parameters are presented for base models (see section 4.2). The unit for v exp is km/sec. 


a change in L* does not affect the shape of the output 
spectra (or color), it only affects the overall energy out¬ 
put. Therefore, we may ignore the effect of the thermal 
pulse on the mean properties (L», T e //) of the central 
stars for a large sample of evolving stars. 

However, a change in the mass-loss rate of the dust 
shell due to the thermal pulse affects the observed SED 
significantly. The various effects of the thermal pulse on a 
dust shell (a superwind and/or chemical transition from 
O to C) in the AGB phase are well investigated (e.g., 
Gronewegen et al. 1995; Justtanont et al. 1996b; Suh 
& Jones 1997). Considering the effects of the thermal 
pulses, we use various mass-loss rates of the dust shell 
(rather than a fixed one) at the end of the AGB phase. 

To make a simplified dust shell model as described in 
section 3.1, we assume that the mass-loss rate and v exp of 
the dust shell at the end of the AGB phase remain con¬ 
stant during the period between the end of AGB phase 
and the start of the post-AGB phase. And we assume 
that AGB dust formation (at T c = 1000 K) continues 
during the period between the end of AGB phase and 
the start of the post-AGB phase because dust formation 
properties during the period can not be easily clarified. 

We assume that the AGB dust formation (at T c = 1000 
K) ceases with the start of the post-AGB phase. t pa 
is defined as the time scale since the start of the post- 
AGB phase. With the start of post-AGB phase, the inner 
radius (Rm) of the dust shell increases with time ( t pa ): 

Rin — Ro A V exp t pa , (5) 

where Rq is R n at t pa = 0. As the post-AGB star evolves, 
T c gets colder as the dust shell detaches. For the dust 
shell, we assume various initial dust optical depths at the 
start of the post-AGB phase so that the corresponding 


mass-loss rates (Equation 3) are reasonably in the range 
of observed values of AGB stars. 

We make various dust shell models with different mass- 
loss rates ( M ) for each core-mass model of the central 
star. For the low core-mass model, we present three 
base models (EM, ML, M) for low-mass O-rich post-AGB 
stars and two base models (CL and C) for C-ricli post- 
AGB stars in Table 3. For the high core-mass model, 
we present three base models (OL, O, and T) for high- 
rnass O-rich post-AGB stars in Table 4. For all dust shell 
models, the expansion velocity ( v exp ) and mass-loss rate 
(M) of the dust shell remain constant during the whole 
post-AGB phase (see section 3.1). 

In Tables 3 and 4, the mass-loss rates are presented for 
base models. We assume that v exp = 15 km/sec for all 
base models. For each base model, dust shell parameters 
(Rim Pin, T io, and T c ) at eight points of time evolution 
in the post-AGB phase are listed. n is defined as the 
initial dust optical depth (no) at the start of the post- 
AGB phase. As the post-AGB phase evolves, no and the 
output parameter T c (the inner shell dust temperature) 
decrease while increases with time ( t pa )■ 

For each base model, we make model calculations for 
various derived models with faster dust shell expansion 
velocities ( v exp = 30, 75, 150, and 300 km/sec). From the 
base model EM, we have two derived models EM2(75) 
and EM2(150) using v exp = 75 and 150 km/sec, respec¬ 
tively. The derived model uses the same model parame¬ 
ters (e.g., Rq and n) and produces the same model SED 
as the base model at t pa = 0. At later post-AGB phases, 
the derived model uses the same stellar parameters as 
the base model but different dust shell parameters (Ri n , 
Pin, and no) because of higher v exp (see Equation 5). 
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Note that M for the derived model (M,j) is given by 

Ma = Mb^, (6) 

15 

where Mb is M for the base model (see Equation 3). 

4.3. Model SEDs 

Figure 5 shows the model SEDs for post-AGB stars 
for six different models. At the beginning of the post- 
AGB phase, the SEDs look like those of AGB stars. As 
the star evolves to a later stage, the SED shows more 
distinct double components, which are typical for post- 
AGB stars. 

Observed SEDs of many post-AGB stars show promi¬ 
nent silicate features (e.g., Hrivnak et al. 1989; Cerrigone 
et al. 2009) as well as the model SEDs for O-rich stars. 
The silicate features become very weak at the end of 
the post-AGB phase especially for the models with high 
v exp . We will compare the theoretical model tracks with 
observations on IR 2CDs in section 5. 

5. COMPARISON ON IR 2CDS 

Figures 1-4 show IRAS , NIR, AKARI, and MSX 
2CDs for AGB stars, post-AGB stars, and PNe compared 
with theoretical models. The small symbols are the ob¬ 
servational data (see section 2.4) and the lines with large 
symbols are theoretical model tracks for AGB stars (sec¬ 
tion 3.1). The lines with filled symbols are theoretical 
model tracks for post-AGB stars (section 4). 

For AGB stars, theoretical models for a range in dust 
shell optical depth is presented (see section 3.1). The 
objects in the upper-right regions on any 2CDs have thick 
dust shells with large optical depths for AGB stars. The 
locations of S stars are restricted to the regions of thin 
dust shells (lower-left regions) on any IR 2CDs. For AGB 
stars, see Suh & Kwon (2011) for a detailed discussion 
on the comparison of the theoretical models with the 
observations on various IR 2CDs. 

To find possible evolutionary tracks from AGB stars to 
PNe on 2CDs, we present the theoretical model tracks 
for the evolving post-AGB phase (see Tables 3 and 4 for 
detailed model parameters). For each post-AGB model 
track, the point at the start of the post-AGB phase is 
indicated by the tip of an arrow and the base model 
name. The base or derived model name is labeled at the 
end point of the model track. 

We find that the theoretical dust shell model tracks 
for post-AGB stars roughly coincide with the densely 
populated observed points of post-AGB stars and PNe 
on various IR 2CDs. Even though some discrepancies 
are inevitable, the end points of the post-AGB model 
tracks are generally converged to the region of the ob¬ 
served points of PNe on most 2CDs. 

If our assumptions on other parameters are right, the 
time scale of the post-AGB phase (832 years) for low 
core-mass models (EM, ML, M, CL, and C) looks rea¬ 
sonable. The models for high core-mass stars (models 
OL, O, and T), for which the time scale is even shorter 
(174 years), match the observations on IR 2CDs fairly 
well with various expansion velocities. 

Unlike AGB stars, C-rich post-AGB stars or PNe typ¬ 
ically show PAH dust emission features at 3.3, 6.2, 7.7, 
8.6, 11.3, and 12.5 /mu (e.g., Justtanont et al. 1996a; Cer¬ 
rigone et al. 2009). This could be due to UV radiation 


from the hot central stars that have evolved past the 
AGB. The observed colors [8.28] —[14.65] and [12] —[25] 
could be bluer than the theoretical model colors which 
do not consider PAH dust. 

The effects of very thin dust shell models (EM, ML, 
and CL) or very thick models (T) can be more clearly 
presented on the IRAS 2CD than any other 2CDs con¬ 
sidered in this work. This is because the [25] —[60] color, 
which uses longer wavelengths, can represent the evolving 
SED better than any other colors (see Figure 5). Other 
model colors can be more easily saturated to certain val¬ 
ues because they can be more severely affected by the 
boundary between two components (the central star and 
the dust shell) and/or silicate dust features. Because 
K s — [12], [9] —[12], and [8.28] —[14.65] model colors are 
rapidly converged to certain values for very small dust 
optical depths, the IR 2CDs in Figures 2-4 are not useful 
to investigate the thin dust shell models compared with 
the IRAS 2CD (Figure 1). Therefore, the model tracks 
for EM2(75), M2(150), ML2(150), ML2(150), ML2(300), 
and CL2(300) are plotted only on the IRAS 2CD (Figure 
1; see Table 3). 

5.1. The IRAS 2CD 

Figure 1 shows the IRAS 2CDs using [25] —[60] ver¬ 
sus [12] —[25] for AGB stars and post-AGB stars com¬ 
pared with theoretical models. The eight regions and 
the curved green line which is the evolutionary track for 
late-type stars were defined by van der Veen & Habing 
(1988). The theoretical model tracks for AGB stars and 
the start points of the post-AGB model tracks roughly 
coincide with the evolutionary track for late-type stars 
which were defined by van der Veen & Habing (1988). 

Brown dashed regions on the IRAS 2CD indicate the 
two groups of post-AGB stars, which are designated as 
LI (Left of IRAS, blue group) and RI (Right of IRAS, 
red group) (Sevenster 2002). Most post-AGB objects are 
distinguished by very red MIR colors (RI region), but 
there is also a group of very blue sources with a strong 
60 fim excess (LI region). A large number of post-AGB 
and PNe samples are mainly distributed in the RI region 
compared with limited samples of the LI region (e.g., 
Yoon et al. 2014). Sevenster (2002) indicated that LI 
objects are more massive and have higher outflow veloc¬ 
ities than RI objects. LI objects may go through an ir¬ 
regular mass-loss phase, preceding the bipolar PN stage. 
On the other hand, RI objects are the traditional post- 
AGB sources (van der Veen & Habing 1988), turned right 
from the evolutionary sequence into IRAS region V and 
beyond. 

The RI objects can be explained by a normal evolu¬ 
tionary sequence for which stars undergo dust shell de¬ 
tachments at the normal AGB termination. For all of 
the post-AGB model tracks with moderate v exp (15 - 75 
km/sec), the end points are located in this RI region. 
This is in accord with the fact that most of the observed 
PNe are converged to the RI region (see Figure 1). This 
trend is more noticeable when we consider only the good 
quality observational data of PNe which are plotted by 
thicker symbols on the 2CD. 

The LI objects could be explained by with extremely 
thin dust shells (models EM2, ML2, and CL2 with faster 
v exp ; see Figure 1). However, the theoretical model used 
for this work may not be applicable because most of the 
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Fig. 5.— Model SEDs for post-AGB stars. We assume that the distance is 1 pc. See Tables 3 and 4 for detailed model parameters. 


LI objects are believed to be non-spherical systems (pre¬ 
cursors of high-mass, bipolar objects or binary systems). 

If we assume that related stars are single and nearly 
spherical objects, the LI objects as well as the objects in 
upper-left region (VI) are likely to be hot objects with 
very thin dust shells. We find that the thin dust shell 
models (EM(2), ML(2), and CL(2)) with various v exp at 
the beginning of the post-AGB phase can explain the 
post-AGB stars in the lower-left region (I, II, HI, and 
VII) on the IRAS 2CD. For the post-AGB stars and 
PNe in the upper-left region (VI and LI), the model re¬ 
quires a very thin dust shell (models EM2, ML2, and 
CL2) with higher v exp (75-300 Km/sec) and a hotter cen¬ 


tral star which is more evolved in the post-AGB phase. 
This requirement of a thin dust shell model may look to 
be in contradiction to the idea that a star at the end 
of the AGB phase would have a very thick dust shell. 
This contradiction could be explained by the effect of 
the chemical transition from O to C in the AGB phase 
(see section 6.3). 

5.2. NIR, AKARI, and MSX 2CDs 

Figure 2 shows IRAS- 2MASS 2CDs using [12] —[25] 
versus K s — [12]. Compared with the IRAS 2CD, the 
boundaries that separate O-rich, C-rich and S stars are 
clearer for AGB stars. The locations of observed PNe 
are converged toward upper-middle region. Unlike other 
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2CDs, the locations of the end points for the post-AGB 
model tracks look to show systematic deviation from the 
densely populated region of the observed PNe. .Kg —[12] 
model colors are generally redder than the observed col¬ 
ors of PNe. 

Post-AGB stars and PNe typically show H 2 line emis¬ 
sion in the K s band (e.g., Davis et al. 2003; Froebrich et 
al. 2011). The 2MASS data in the I\ s band (2.17 /jrn) 
may show this effect. If the effect of H 2 emission on the 
K s band is stronger than the one from the PAH emission 
on the [12] band, the observed K s — [12] colors could be 
bluer compared with the theoretical model colors which 
do not consider the H 2 line emission. Another possible 
reason for the difference could be that the 2MASS flux 
may include only the central star while the IRAS fluxes 
may include the extended post-AGB dust envelope (see 
section 2). 

Figure 3 shows AKARI-IRAS 2CDs using [9]-[12] 
versus [12] —[25]. Generally, the end points of the O-rich 
and C-rich post-AGB model tracks are converged to the 
region of observed PNe. Compared with other 2CDs, 
it is difficult find any boundaries that separate O-rich 
AGB stars and C-rich AGB stars. But the boundary 
that separates AGB stars and PNe is clear with the aid 
of the distinguishing [12] —[25] color. The model [9] —[12] 
colors can be affected by some dust species (alumina, 
SiC, and PAH) which are not considered in this work. 

Figure 4 plots MSX 2CDs using [8.28]-[14.65] ver¬ 
sus [14.65] —[21.34]. Unlike other 2CDs, the base mod¬ 
els (OL, O, and T) for high core-mass stars can match 
the observations better than the models with faster v exp . 
Generally, this MSX 2CD looks to require a faster evo¬ 
lution of the post-AGB phase than other 2CDs. The 
[14.65] —[21.34] color is useful to distinguish different 
classes of objects on the 2CD. For the [8.28] —[14.65] 
color, the PAH emission feature at 8.6 /rrn from post- 
AGB stars and PNe could make the observed colors bluer 
(downward on the 2CD) compared with the theoretical 
models. 

6. DISCUSSION 

6.1. Uncertainties and alternative models 

To investigate spectral evolution from AGB stars to 
PNe, we have performed radiative transfer model calcu¬ 
lations for the evolving dust shells and central stars in 
the post-AGB phase. We have assumed that the central 
star emits blackbody radiation. Uncertainties remain in 
a number of aspects. 

The spherically symmetric dust shell models used for 
this work could be too simple for post-AGB stars because 
their shapes are often believed to be non-spherical (e.g., 
bipolar or disk-like objects). However, the theoretical 
models for non-spherical dust envelopes (disks or toruses) 
which are applicable to the highly non-spherical objects 
should have a much larger number of assumed model 
parameters which can be difficult to be compared with a 
large sample of observed objects. 

Our assumption of the blackbody radiation for the cen¬ 
tral star inside a dust shell could be too simple espe¬ 
cially for post-AGB stars for which gas-phase radiation 
from the outer shell can be important. For the evolu¬ 
tion of a post-AGB star with an expanding gas shell, 
Volk (1992) used the output of the photoionization code 


CLOUDY (Ferland 1993) as input for a dust shell model 
(CSDUST3). But this could have a problem of incon¬ 
sistent treatment of dust temperatures, van Hoof et al. 
(1997) used the CLOUDY code considering dust emis¬ 
sion more consistently for an optically thin dust shell. 
However, most dust shell models at the beginning of the 
post-AGB phase require large dust optical depths. 

For dust opacity, we have not considered some dust 
species: alumina, SiC, and PAH. Suh & Kwon (2011) 
used a mixture of alumina and silicate for O-rich AGB 
stars and a mixture of SiC and AMC in modelling C- 
ricli AGB stars for various model tracks. We also used 
the same mixtures for testing post-AGB models, which 
produced different model tracks but the overall match 
to the observations looked similar. Unlike AGB stars, 
C-rich post-AGB stars or PNe typically show PAH dust 
features. We have discussed possible effects of PAH on 
IR 2CDs in section 5. 

6.2. Comparison with other works on post-AGB stars 

Even though both our dust shell model and van Hoof 
et al. (1997)’s work assume a spherically symmetric con¬ 
tinuous dust shell, van Hoof et al. (1997) used signifi¬ 
cantly different schemes in their radiative transfer code 
(CLOUDY) from those used for this work (RADMC-3D). 
The strong point is that their model can consider some 
gas-phase radiation (bound-free emission for hot mod¬ 
els). As commented in the paper, the weak point is that 
the CLOUDY model for a thick dust shell could be unre¬ 
liable. van Hoof et al. (1997) provided many meaningful 
results, but some of their results need to be considered 
more carefully because dust optical depths at the begin¬ 
ning of the post-AGB phase are generally large. 

On the IRAS 2CD, van Hoof et al. (1997) presented 
the red loop model track which causes the [12]-[25] color 
to loop back towards the blue, before the star resumes 
its normal evolution to the red. Some of our models (O 
and T) also show similar tracks, but they do not cover 
deep into the LI region. Sevenster (2002) tried to explain 
the LI objects on the IRAS 2CD (see section 5.1) using 
the leftward model track in the early phases of the red 
loop track. However, their dust shell model (CLOUDY) 
could be unreliable because the corresponding dust opti¬ 
cal depths (for the leftward model track) are very large. 
Volk (1992) and Ortiz et al. (2005), who used the out¬ 
put of the photoionization code (CLOUDY) as input for 
a dust shell model (CSDUST3), also did not find the red 
loop model track. 

If our assumptions are right, the LI region of the IRAS 
2CD can be covered by the model tracks with thin dust 
shells at the starts of their post-AGB phases and larger 
detachments which can be achieved by higher expan¬ 
sion velocities (e.g., models EM2, ML2, and CL2 in this 
work). These models make upward model tracks rather 
than leftward (see Figure 1). Of course, any spherically 
symmetric dust shell models would be inappropriate for 
highly non-spherical objects in the LI region. 

For the low core-mass model, we also used the model 
for M core = 0.605 Mq with the post-AGB mass loss rate 
at the standard value (see section 4.1). The correspond¬ 
ing dynamical time scale of the post-AGB phase is 2093 
years. This model produced the model results which are 
too far off from the observed points of post-AGB stars 
and PNe on all 2CDs. Volk (1992) and Ortiz et al. (2005) 
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investigated spectral evolution of post-AGB stars using 
the similar scheme for the evolution of the central star. 
They also showed that a fast evolution (the time scale 
of less than about two thousand years) was required to 
match the observations on 2CDs. Because they assumed 
an earlier start of the post-AGB phase, the dynamical 
time scales from Ortiz et al. (2005) are generally larger 
than those used for this work. 

6.3. The chemical transition in the AGB phase 

As stars evolve into the thermal pulsing AGB phase, 
the abundances of some elements in the stellar atmo¬ 
sphere may change by the episodic third dredge-up pro¬ 
cess after each thermal pulse. When AGB stars of in¬ 
termediate mass range go through carbon ’dredge-up’ 
processes, and thus the abundance of carbon is larger 
than that of oxygen, O-rich dust grain formation ceases 
and the stars become visual carbon stars. After that 
phase, carbon-rich dust grains start forming and the stars 
evolve into infrared carbon stars with thick C-rich dust 
envelopes and very high mass-loss rates (e.g., Iben 1981; 
Chan & Kwok 1990; Suh 2000). 

S stars are generally regarded as intermediate between 
M-type and carbon stars in their properties. However, 
this M-S-C evolutionary sequence can be different de¬ 
pending on the mass and metallicity; some stars may 
remain in S-type until a next thermal pulse but some 
stars may skip the S-type star phase (e.g., Gronewegen 
et al. 1995). SC or CS stars are likely to be in a transi¬ 
tional phase from an S star to a carbon star. Generally, 
dust formation in envelopes of S stars is less efficient be¬ 
cause of the lack of free O or C to form dust. Most S 
stars show very weak dust (silicate or carbon) emission 
features and low mass-loss rates (e.g., Hony et al. 2009; 
Smolders et al. 2012) and their locations are restricted 
to the regions of thin dust shells (lower-left regions) on 
any IR 2CDs (see Figures 1-4). 

To reproduce the post-AGB stars in lower-left and 
upper-left regions on the IRAS 2CD, we need to use 
very thin dust shell models with higher v exp at the be¬ 
ginning of the post-AGB phase (see section 5.1). S stars 
or visual carbon stars which have very thin dust shells 
are believed to be products of the chemical transition 
from O to C in the AGB phase. The objects which be¬ 
come S stars or visual carbon stars nearly at the end of 
the AGB phase are likely to have very thin dust (silicate 
or AMC) shells at the beginning of the post-AGB phase. 
These objects with with higher v exp could be the post- 
AGB stars or PNe in the upper-left region (VI and LI) 
on the IRAS 2CD. 

7. SUMMARY 

We have presented various IR 2CDs for AGB stars, 
post-AGB stars, and PNe using catalogs from the avail¬ 
able literature for the sample of 4903 AGB stars (3373 
O-AGB; 1168 C-AGB; 362 S-type), 660 post-AGB stars 
(326 post-AGB; 334 pre-PNe), and 1510 PNe in our 
Galaxy. For each object in the catalog, we cross-identify 
the IRAS, AKARI , MSX, and 2MASS counterparts. 
For the large sample of stars, we have presented various 
IR 2CDs using the IRAS PSC, AKARI PSC, MSX 
PSC, and NIR (2MASS data at K s band) data. 


To find possible evolutionary tracks from AGB stars 
to PNe on the 2CDs, we have made simple but reason¬ 
able assumptions on the evolution of the dust shell and 
the central star in the post-AGB phase. We assume that 
AGB dust formation (at T c = 1000 K) ceases when the 
post-AGB phase starts. The dust shell detaches with the 
start of post-AGB phase because there is no more dust 
formation. We have performed radiative transfer model 
calculation for detached dust shells around evolving cen¬ 
tral stars in the post-AGB phase. 

We have found that the theoretical dust shell model 
tracks using dust opacity functions of amorphous sili¬ 
cate and amorphous carbon roughly coincide with the 
densely populated observed points of AGB stars, post- 
AGB stars, and PNe on various IR 2CDs. Even though 
some discrepancies are inevitable, we have found that 
the end points of the theoretical post-AGB model tracks 
are generally converged to the region of observed points 
of PNe on most 2CDs. The discrepancies could be due 
to the limitation of the theoretical model used for this 
work. The dust shell model did not consider some dust 
species (e.g., PAH) and gas-phase radiation processes. 
The spherically symmetric model would not be applica¬ 
ble to a major portion of the observed objects which are 
highly non-spherical. 

In this work, comparison of the model results with the 
observations on the IR 2CDs have indicated that the du¬ 
ration of the theoretical post-AGB phase needs to be 
rather short (within a thousand years) to match the ob¬ 
servations. This could be because relatively small mass 
stars are not observed as post-AGB stars or PNe. The 
duration of the post-AGB phase for very low-mass stars 
can be as long as 10 4 -10 5 years (see section 4.1). A major 
portion of these very low-mass stars may evolve to the 
white dwarf stage without ever becoming a PN. 

We have discussed two sequences of the post-AGB evo¬ 
lution on the IRAS 2CD, which are designated as LI 
(blue group) and RI (red group). The red objects are 
likely to undergo dust shell detachments at the normal 
AGB termination. Most of the observed PNe are con¬ 
verged to the RI region as well as the results of all the¬ 
oretical dust shell models with moderate v exp . On the 
other hand, most of the blue objects in the LI region are 
likely to be non-spherical systems (precursors of bipo¬ 
lar PNe or binary systems) to which the theoretical dust 
shell models may not be applicable. However, we have 
demonstrated that it is also possible that the blue post- 
AGB stars in the LI region as well as in the upper-left 
region of the IRAS 2CD have very thin dust shells with 
higher v exp because they had thin dust shells at the starts 
of their post-AGB phases. We have argued that an ob¬ 
jects in its early post-AGB phase may have a very thin 
dust shells if the chemical transition from O to C due to 
a thermal pulse occurred in its late AGB phase. 
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